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Abstract. In order to explain the recently reported peak at 3.55 keV in the galactic x- 
ray spectrum, we propose a simple model. In this model, the Standard Model is extended 
by including a neutral spin-3/2 vector-like fermion that transforms like a singlet under SM 
gauge group. This 7.1 keV spin-3/2 fermion is considered to comprise a portion of the observed 
dark matter. Its decay into a neutrino and a photon with decay life commensurate with the 
observed data, fits the relic dark matter density and obeys the astrophysical constraints from 
the supernova cooling. 
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1 Introduction 

Recently X-Ray emission at ~ 3.55 keV has been observed in the XMM- Newton X-Ray 
observatory [1, 2] in many Galaxy clusters and in the Andromeda Galaxy spectra. The 
observed flux and any X-Ray line energy measured in the MOS spectra is given by 

^MOS ^ X 10"® photons cm'^ sec"^ 

^ 3 54 ^ Q_Q2 keV (1.1) 

The source of this line is yet to be identified. An attractive possibility, considered in the 
literature to explain the observed flux and energy, is to attribute it to the decay/ annihilation 
of some dark matter particle which is stable over cosmological time scale and can account 
for at least a significant fraction of dark matter relic density with mass and decay life time 
commensurate with the observed data. Sterile Neutrino of mass 7.1 keV capable of produc¬ 
ing warm Dark Matter (WDM) density through resonant or non-resonant production with 
parameters required to produce the observed signal is an attractive proposition discussed in 
the literature [3-7]. R-parity violating decays of the lightest super-symmetric particle (LSP), 
the decay of gravitons and axions, into neutrino photon pair and the decay of scaler field cj) or 
axion-like pseudoscalar fields a into photon pairs as possible explanation of the signals have 
been considered in the literature [7-16] with varying success. The scalar case is of particular 
interest because the scale of the new physics may involve super-symmetry (SUSY) which 
conforms to the expectations from the physics of moduli. 

Several new physics models beyond the standard model (SM) predict the existence of 
spin-3/2 particles. In models of super-gravity, the graviton is accompanied by spin-3/2 grav- 
itino super partner. In models of composites [17], the top quark has an associated spin-3/2 
resonance. New physics models may include exotic fermions and gauge bosons which are 
not present in the SM. Spin-3/2 fermions also exist as Kaluza -Klien modes in string theory 
[18, 19] if one or more of compactification radii are of the scale lower than the Planck scale. 

The prediction of spin-3/2 particle as a cold dark matter has been made by several 
authors in SUSY models [20, 21]. Gravitinos with mass in the keV range have been studied 
as the probable WDM candidate in various SUSY models [22-25] even before the observation 
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of 3.55 keV X ray emission. Recently, authors of the reference [26, 27] have studied the 
implication of the effective four fermion interactions involving the DM spin-3/2 particle on 
relic density, the antiproton to proton flux ratio in cosmic rays, and the elastic scattering 
off nuclei ( direct detection) in the effective field theory approach. Constraints from direct 
detection of dark matter exist in literature on spin-3/2 WIMP candidates [28]. 

A recent comprehensive analysis by the authors of reference [29] demonstrated that 
the measured flux of the 3.55 keV line can be accounted for, by the conventionally known 
plasma lines without invoking the dark matter decay as its origin. This explanation, however, 
requires the fixing of the abundances of different elements which are still uncertain to a certain 
degree. We, thus, feel that it is worthwhile to investigate alternative interpretations that are 
consistent with the other astrophysical and cosmological data. 

In this paper, we consider a new neutral spin-3/2 fermion assumed to be a vector-like 
SM singlet. We will consider the decay of this 7.1 keV DM particles into a neutrino-photon 
pair (x with decay life commensurate with the observed galactic X-ray spectrum. 

This spin- 3/2 particle could exist as fundamental particle or could be a bound state of SM 
neutrino and U{1) gauge bosons. We will explore the possibility of such an exotic spin-3/2 
particle to constitute the relic dark matter for a reasonable choice of parameters and confront 
the model from cosmological and astrophysical constraints. 

In section 2, we describe the spin-3/2 fermion model. In section 3, we discuss the 
implication of the model to explain the observed galactic X Ray spectrum data. In section 
4, we obtain the relic abundance and the resulting constraints on the model parameters. In 
section 5, we discuss the bounds obtained from supernova energy loss. Section 6 is devoted 
to results and discussion. 


2 The spin-3/2 Model 


The standard model is extended by including a spin-3/2, vector like particle Xi whose right- 
handed (RH) as well as left-handed (LH) projections transform the same way under SU{2) x 
SU{1). We further let x to be a SM singlet. Spin-3/2 free Lagrangian is given by 


^ = Xim{Px) Xu{Px) where 

= {i ^ — i ( 7 ^ 7 ^ -|- 7 ^ 7 ^) + *7^ ^ 7 ^ + 

Here Xfi satisfies A^'^Xu = 0. For on mass-shell X: we have 

X^XfiiPx) = 0 = d^XfiiPx) = mx)XfiiPx)- 

The spin-sum for spin-3/2 fermions 

3/2 3/2 

•S^f^,^{p) = u'f,{p) Kip) and S ~^ v^p) Kip) 

i=-3/2 i=-3/2 


are given by 


S%{p) = rux) 


1 2 1 , / 

Pfiu - 77m lu - j^PfiPu =F 7 : 7 ^ Kt^Pu - y^Pti) 


X 


3m. 


( 2 . 1 ) 


( 2 . 2 ) 


(2.3) 


(2.4) 
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respectively. 






The most general leading order standard model gange invariant interaction between 
the spin-3/2 SM singlet x and SM spin-1/2 fermions is given by the effective dimension six 
operators: 


/-Eff. 


A2 

2=1 




^ A 7 " [ 7 ^, l'"] + § 4 7 “ x'^d) B 


A2 


A2 




(2.5) 


where = idf, - i {gs/2) A“G“ - i {g/2) - ig'Y = {cosOF^^ - sinOZ^^), 

(j) = iT 2 (j) and is the SM lepton donblet. The weak 17(1) hyper-charge Y for cj) and x are 
1/2 and 0, respectively. 

in view of the on-mass shell conditions as given in Eq. (2.2), the second operator O 2 
vanishes and the third operator O 3 becomes identical to the first operator Oi. Therefore we 
are left with only one coupling constant C, which can be simplified to give after symmetry 
breaking: 


Bi 


Cvo 

A2 


l^kL {PuJl/iXp (Px) (cos 0 


smez^^p). 


Here vq is the SM Higgs vacnnm expectation valne and A is the new cnt-off scale. 


( 2 . 6 ) 


3 Galactic X-ray spectrum 


The decay width for x 


t'e 7 is given by 


Cvq cos 6 
A2 


- 

IGvr ^ 


27 

■ C ■ 

2 

■ A ■ 

-4 

[ "ix ] 


10-9 


100 TeV 


WkeVJ 


sec 


-1 


(3.1) 


Here we have taken the conpling of spin-3/2 particle with only one generation (say for the 
first generation only) of SM nentrino. 

The expected X-ray flnx is proportional to the density of the decaying dark matter x- 
The WDM which in the case considered here constitntes of spin-3/2 SM singlet, is believed to 
comprise a portion of the observed DM relic abnndance with CDM as the dominant component 
[30, 31]. If the X-ray galactic signal is interpreted as coming from the spin-3/2 WDM x 
decaying into a nentrino and a photon pair, the reqnired valne of the life time of x shonld be 
given by 1.4 / X 10^® seconds, where / (0 < / < 1) is the fraction of the relic dark matter 
density contribnted by the WDM x- At / = 1.0 the WDM x wonld acconnt for the entire dark 
matter relic density with a choice of new physics scale A of the order of ~ 100 TeV along with 
the conpling constant C ~ 10“®. The small valne of C shonld not be snrprising as it can be 
considered to be a measnre of trilinear lepton nnmber violating conpling and hence natnrally 
small. Similar sitnation occnrs in snper-symmetric models of R-parity violating interactions 
considered in the literatnre [9-16] as possible explanation of the observed galactic X-ray flnx. 
In realistic model, the mixing between photino and nentrino for example, is snppressed by a 
small parameter ~ 10“^*^ characterising lepton nnmber violation [14, 32-34]. 
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If, the 7.1 keV signal, on the other hand, is interpreted as coming from pair annihilation 
of 3.55 keV spin-3/2 DM into two photons, the annihilation cross-section {av)ann. has to 
match with the best-fit decay-width of 7.1 keV DM i.e. 

{<rv).nn. ~ (3.2) 

nx 

where = p^jm^ ~ (lO^ — 10^) cm“^ is the number density of spin-3/2 DM. This trans¬ 
lates into (cru)ann.fit — 2 X 10“^® GeV“^. 

The spin-3/2 particles can couple to two photons through 17(1) gauge invariant dimension 
seven effective Lagrangian 

Ant. = (3.3) 

This gives an annihilation cross-section cr(xx 77) ~ m^/(7rA®) and the desired annihi¬ 

lation rate is achieved for A < 0(100) MeV (for < 1), which is clearly unphysical. Thus, 
it is unlikely that the observed galactic X-ray signal can be explained by DM y’s annihilation 
into photons. 

4 Relic Abundance of Spin-3/2 

Since the x’s couple weakly to the SM particles and are nearly stable with a lifetime com¬ 
parable to the age of the Universe if they have to account for the observed X-ray flux, they 
will decouple early when they are relativistic. They will, therefore, contribute to the present 
mass density of the Universe as DM. Their abundance at decoupling is nearly equal to the 
photon density at that time. During the adiabatic expansion of the Universe, their num¬ 
ber densities remain comparable. A rough estimate of the bound on y mass can be ob¬ 
tained just like the bound on the neutrino mass [35] by requiring that the ratio of DM x 
density to the critical density remain less than one. This gives < 12.8g* (Tu) /g^ff 
eV. The effective number of degrees of freedom g* (Td) at decoupling time of electroweak 
symmetry breaking transition is found to be about 113.75. In the computation of g* {Td), 
we have included the effective degrees of freedom from all SM particles and y, x spin-3/2 
DM particles. However, in the MSSM, g* [To) is much larger ~ 228.75 and thus it is not 
reasonable for the spin-3/2 DM particle y to have a mass of the order of about 7.1 keV. 

The relic abundance of dark matter x de¬ 
pends on the sources of production of x 
in the early Universe. The leading order 
processes (shown in Figure 1) that main¬ 
tain the DM X ill equilibrium with the rest 
of the SM plasma are the decay rate of 
Z —>■ xT'e + X^^e and the 2 —>■ 2 pair anni¬ 
hilation rates, namely, V{Z —?• xA + X^e)) 

XA + X^e} and ct(V+1U" 

XT'eFX^e) where Sj. means summation over 
all SM fermions (quarks and leptons). Us¬ 
ing the interaction Lagrangian given in Eq. 
(2.6), the decay and spin averaged annihi¬ 
lation cross-sections can be computed in a 


(i>.) 




{^) i>e (i^e) \y+ 


x(x) 




i>e (Ve) 


Figure 1. Relevant diagrams for decay and •pro¬ 
duction of the DM candidate x- 
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straightforward manner. We obtain: 
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a{W+W- 


XV e + XVe) 


A'KaC'^VQ 1 
A^ 2887r 


1 - 
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mr 
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s (s — m |)2 + r^m^ 
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I 18 m| 
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TTIjj/ 772x1/ 

1 + 20—54^ + 12^ 


(4.3) 


where gy and 5 :^^ are the vector and axial vector couplings of respective fermions in SM. 

If the decay and annihilation rates are much smaller than the Hubble expansion rate at 
the temperature of the order of Elctro-Weak (EW) symmetry-breaking scale, the spin-3/2 DM 
particle x will never be in thermal equilibrium. The Hubble expansion rate at a temperature 
T is given by H{t) ~ 1 . 665 *r^/mpi where g* is the effective number of relativistic degrees of 
freedom at the temperature T. The decay Z -A- xVe + comes into play only below EW 
symmetry breaking phase transition temperature Tew ~ 150 GeV. The Z bosons go out of 
the equilibrium below roughly 5 GeV, the other SM fermions remain in equilibrium much 
below this temperature. 

The decay and annihilation rates can be estimated from Eqs. (4.1)-(4.3). The Z decay 
rate is given by 


^{Z -A XVe + XVe) 


7.7 X 10^^ 


A 


1 -4 


IGeV 


GeV. 


(4.4) 


The leading terms in the cross-section corresponding to x production through // annihilation 
and Wfusion processes are given by 


O-ififi XVe + XVe) 


2.3 X 10^^ 


s 


■ A ■ 

[iGeV^J 


IGeV 


GeV-2 
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and 


a{W~^W —> x’^e + X^e) ~ 6.96 x 10^ 
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3 

■ A ■ 

llGeV2j 


IGeV 


1 -4 


GeV 


-2 


(4.6) 


respectively. One can obtain the constraint on the effective coupling C by demanding the 
thermal average r(Z x^^e+X^'e), (^cr (j2f fl ^ X^e + and {a{W^W~ x^^e + Xt^e)\nv) 

to be less than the H(T) for T ~ 150 GeV (ie. at the EW phase transition temperature). 
Therefore, using g* = 113.45, we obtain 


C 


and 


C 


C 


A 


IGeV 

A 

IGeV 


A 

IGeV 

-2 


1 -2 


-2 


: 10 from Z decay. 

(4.7) 

from 

(E/// X^e + X^e) , 

(4.8) 

from 

W^W~ -z Xi^e + X^e- 

(4.9) 


The thermal averaged cross-sections (cr|nf) are estimated using the relation s = 4 {E) where 

(E) = 3.15 T and 2.7 T, and n = - and ^^^gT^ for Fermi-Dirac and Bose-Einstein 

particles, respectively. 

The relic density of the spin-3/2 DM x can be evaluated by solving the Boltzmann 
equation for the evolution of the number density of DM x and is given by 


ri^ -t- = -{T{Z xz^e))(n^ - Uq) - E/ ( cr | ^ // -^ XJ"e ) \v) - n^^) 


— (cr(iy'“lT xt^e)\v) {n 


X^Ue 




(4.10) 


Here, Uq is the equilibrium number density of species i. The region of validity of the equation 
is when all the SM particles are in thermal equilibrium unlike the DM candidate x which 
is realized for T < (15 — 20) rrii. We can than put = 0 in the R.H.S. of this equation. 
Ghanging the variable from time to temperature, the equation can be put in the form: 


djx 

dz 


r(-Z' —)• xw) sr^ X^e)\v) I fi -^2 

^Jo ' / ^ KZ^ ’ 


Km?^ 


(4.11) 


where z = mz/T, = n^/T^, /* = Hq/T^, and K = /mpi. We use Boltzmann 

distribution functions for both the fermions and bosons, i.e. 

fiirm/T) = /• (^^) = ^ J"p\-z/P'+‘P'dp. (4.12) 


The thermal averaged decay rate and annihilation cross-sections can be expressed, following 
Ref. [36, 37], as 

( ‘z') 

(4.13) 


(r(z —x^e)) — r(z — x^e)- ^ ^ 


K2 {z) 
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and 


W.n ^ Xt'e)) = (y) ''(/./i ^ Xi'Hs - 4mf)v/iA'. *. (444) 

Here, Ki^ 2 ix) are the modified Bessel functions. In terms of scaled number density defined 
as = f^Kmz and by using the expressions for thermal averaged decay width and the 
annihilation cross-sections given in Eqs. (4.13) and (4.14), the Boltzmann equation can be 
written as 


dz 


= V{Z 




1 


K 2 {z) mz 


zfSiz) 


1 z 




8 xf K^ixiz) 


{y - ^xl)y/yKi{y/yz)a{fifi 


X^e)^ {fo{xiz)f 


dy. 


(4.15) 


We solve the above Boltzmann equation for the scaled number density of spin-3/2 dark 
matter particle x for 0.6 < z < 18 corresponding to 5 GeV < T < 150 GeV. 

The contribution of the proposed 7.1 keV spin-3/2 fermion x to the relic dark matter 
density is obtained by numerically solving the Boltzmann Eq. (4.15) from the electroweak 
phase transition temperature to the freeze-out temperature of Ws and Z's. The scaled 
number density N^'s for the leading processes that maintain the dark matter x in equilibrium 
with the rest of SM plasma are obtained to be 


Ny, (r(z)) ~ X 10 


19 


A 


IGeV 




X 10 


20 


A 


IGeV 


-4 


and {a{W+W-)) ~ x 10 


23 


A 


1 -4 


IGeV 


(4.16) 


for the Z-decay, fermion-antifermion annihilation and fusion processes, respectively. 

We hnd that the contribution of spin-3/2 DM fermion to the relic density from 
boson fusion process is about three order of magnitudes greater than the contribution from 
the rest of the processes. For our estimate of the dark matter density, we use the N^. ~ 
X 10^^ [A/ 1 GeV]~^. Thus, the number density of x nt the electroweak phase transition 
temperature is given by n^/ ~ V^/ (AT M^) ~ 10^^/ {K M^). The number density 

of x’s as the Universe cools to the present day is estimated to be ^xIto Ti X Id^^liCKmz) 
where Tq is the present day temperature (Tq = 2.7377) and = g*{TEw)/9* (To) ~ 33.85. 
The present day dark matter relic density PxItq ~ x 10“® [A/lGeV]~^ GeV^ is then 
obtained by multiplying the number density with its mass m-^. Since, the critical dark 
matter density rsj § ^ h? 10 GeV^, the = Px!Pxc is computed as 




0.11 X 10“^^ 


A 

IGeV 



(4.17) 


However, the desired value of ~ ^dm = 0.11 will be obtained for C [A/IGeV] ^ 

10-21. 
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5 Supernova Energy Loss 


The 7.1 keV spin-3/2 dark matter y can be a source of significant energy loss in the supernova 
core. The emission rates for SN 1987 A have been extensively studied for weakly interacting 
DM candidate particles like axions, gravitinos, right handed neutrinos, majorons, low mass 
neutralinos, Goldstone bosons etc. in new physics models. Constraints have been put on the 
properties and interactions of these particles [38-44]. The SN bound on neutrino magnetic 
dipole moment have been one of the tightest [45]. In our estimate of the constraints on the 
parameters of our model, we would use the Raffelt criterion [46] that new source of cooling 
should not exceed the emissivity e/p = 10^® ergs per gm per sec. The main source of y pair 
production in the core of SN is through the yug and/ or yug production processes. The 
emissivity i.e. the energy emitted per unit time and volume, is 



e = 


where |M|^ is the matrix element squared, summed over the initial and final states and 
fi = [ex.p{Ei — iJ,i)/T + 1]~^ is the Fermi-Dirac distribution for the ith particle. 


In the supernova core immediately after the collapse, the temperature is high being of 


the order of tens of MeV. Even though the nucleons are nearly non-degenerate, the electrons 


are degenerate and the neutrinos are trapped. The core has a fixed value of the lepton number. 
Thus, there also exists a sub-dominant energy loss process via the neutrino pair annihilation 
yz 4 -|-yug. Since, the coupling of the dark matter particle y to SM fermions is extremely 
weak, the y’s once produced freely stream out of the SN core, their mean free path being 


greater than the core radius. We thus have ~ 0. Carrying out the phase space integrals 


and making a change in variables from Ei, E 2 , 9 to E^ = Ei + E 2 , E- = Ei — E 2 and 
s = 2mg -|- 2 E 1 E 2 — 2 pi.p 2 cos 6, we get 



In deriving the above expression, we have neglected the Pauli blocking terms for the 
final state particles y and Ug which is an excellent approximation for Ug, y and y. We have 
similar expression for the process UgUg —>■ yug-l-yug. The cross-section for these processes has 
been evaluated in Eq. (4.2). 


The core density lies anywhere between 3 x 10^^ to 3 x lO^^gm/cc. At a core temperature 


of about 40 MeV, the electron chemical potential is //g ~ 200 MeV and |Ug — ~ 50 MeV. 


In our estimate of the energy loss, we consider the core density to be 3 x 10^“^ CeV with a 


core temperatures 30 (50) MeV and electron and neutrino chemical potentials 200 (150) MeV 
and 150 (100) MeV, respectively, and evaluate the energy loss integral numerically. 

Constraints from supernova cooling are obtained by numerical integration of the emissiv¬ 


ity expression (5.2) for the process e"*“e —)> yug -|- yug at T = 30 MeV and electron chemical 
potential /Xg = 200 MeV, we obtain 


e(e+e yug -F yug) 


= 2.2 X 10®^ 


- 8 


Pcore 


IGeV 


ergs/gm/s 


(5.3) 












where we have taken the core density pcore to be about 3 x 10^'^gm/cc. Requirement of 
< 10^®ergs/gm/cc constrains 


pcore 


C 


A 


IGeV 


-2 


< 6.7 X 10“^®. 


(5.4) 


Core temperature of 50 MeV and = 250MeV results in a somewhat tighter constraint 
C [A/ 1 GeV]~^ < 10“^^. The contribution from the process t'gTe —^ X^e + X^e is totally 
negligible being roughly 10 orders of magnitude smaller compared to the annihilation process. 


6 Results and Discussion 

6.1 Summary 

We summarize the constraints on the parameters of our 7.1 keV spin-3/2 dark matter particle 
y from cosmological and astrophysical observations. We observe that the ratio of the coupling 
and the square of the cut-off scale Cjh? associated with spin-3/2 particle y of mass 7.1 keV 
can be constrained as 

• C Qev ] ^ ^ 10“^^ from the consideration of y as a WDM candidate accounting for 

the entire observed relic dark matter density = ^dm = 0.11, 

• C [ iGev ] ^ ^ ^ 10“^® from the the rapid cooling of the supernova through the 

emission of y and, 

• C [ iQev ] ^ ~ ^ 10“^^ from the lifetime of y through its decay y —)• VeX- 

These combined constraints on the parameter space of coupling C and the cut-off scale A 
arising from its appropriate lifetime, contribution to relic density and supernova cooling are 
shown in hgure 2. The curves marked / = 0.01 and / = 1 correspond to the life time 
required for the observed X-ray flux for WDM y contribution = 0.01 x VLum and 

= ^dm respectively. We hnd that the constraints from the cooling of supernova 1987A 
and the DM relic density = 0.11 enclose an allowed band (shaded with yellow lines in 

the hgure) in the parameter space spanned by C and A. The parameter space shaded in green 
is forbidden. 

We thus see that a minimal extension of the SM by adding a spin-3/2 SM singlet with 
mass 7.1 keV can account for the dark matter in the Universe, while at the same time ex¬ 
plaining the 3.55 keV X-ray line in the galactic X-ray spectrum through its decay y —>■ 1/7. 

6.2 Outlook 

Recently, superconducting detectors are proposed for direct detection of light DM particles of 
mass as low as 1 keV through electron recoil from DM-electron scattering in superconductors 
[47]. It will be worthwhile to study the DM model discussed in this article to compute the DM 
scattering rates with electrons in a superconducting environment where electrons are highly 
degenerate and the scattering is inhibited by the Pauli blocking and to explore the feasibility 
of detecting the proposed DM particle. We leave this for the future work. 
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Figure 2. Combined constraints on the coupling C and cut-off scale A from contribution to the 
relic density as WDM, the rapid cooling of supernova SN 1987 A and decay of the spin-3/2 particle. 
Allowed region (with yellow lines) of the parameter space is bounded by the maximum DM relic density 
constraint = 0.11 and from supernova cooling of SN 1987 A. The region in green is forbidden. 
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